Model for the Scaling of Stresses and Fluctuations in Flows near Jamming 
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We probe flows of soft, viscous splieres near tfie jamming point, which acts as a critical point 
for static soft spheres. Starting from energy considerations, we find nontrivial scahng of velocity 
fluctuations with strain rate. Combining this scaling with insights from jamming, we arrive at an 
analytical model that predicts four distinct regimes of flow, each characterized by rational-valued 
scaling exponents. Both the number of regimes and values of the exponents depart from prior 
results. We validate predictions of the model with simulations. 

PACS numbers: 47.57.Bc, 83.50.Rp, 83.80.Iz 



The past few years have seen enormous progress to- 
wards understanding the static, "jammed" state that oc- 
curs when soft athermal particles are packed sufficiently 
densely that they attain a finite rigidity [l -Sj. Such sys- 
tems may flow when shear stresses are applied, and in 
seminal work, Olsson and Teitel addressed the relation 
between strain rate, shear stress and packing fraction in 
a simplified numerical model for the flow of soft viscous 
spheres [4]. When rescaled appropriately, the data for 
strain rate 7, shear stress a and packing fraction (jj were 
found to collapse to two curves, reminiscent of second 
order-like scaling functions, and a large length scale was 
found to emerge near jamming. Since then, qualitatively 
similar results have been obtained in simulations of a 
number of flowing systems [5-1)] , but there is little agree- 
ment on the actual value of scaling exponents, nor on the 
relation to jamming in static systems. 

Here we describe an analytical model that connects the 
scaling of static systems to the scaling of both the velocity 
fluctuations and the shear stress of flowing systems near 
jamming. The model is built around a "viscoplastic" ef- 
fective strain 7og = + 7(jy„, where 7dyn is a dynamic 
contribution set by the strain rate, and 7y stems from the 
(dynamical) yield stress and is controlled by the distance 
to jamming. We show that steady state power balance 
dictates nontrivial scaling of 7dyn with strain rate, and 
propose a nonlinear stress-strain relation that leads to a 
closed set of equations predicting a rich scaling scenario 
for flows near jamming. We verify central ingredients of 
the model and our predictions for the rheology numeri- 
cally in Durian's bubble model for foams |TT] . Our simple 
model captures and predicts the rheology and fluctua- 
tions starting from the microscopic interactions; it also 
indicates the need for, and provides, new ways to present 
and analyze rheological data near jamming. 

Numerical Model — The two-dimensional Durian bub- 
ble model stipulates overdamped dynamics in which the 
sums of elastic and dissipative forces on each bubble, 
represented by a disk, balance at all times [H]. Forces 
are pairwise and occur only between contacting bubbles. 
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FIG. 1: Viscous stress (a) and elastic stress (b) for packing 
fractions = 0.80, 0.8424, and 0.87. The dashed line ~ 7° ''*. 
CTvisc scales linearly with 7 (solid line) and lies below the elastic 
stresses for 7 < 0(10"^). 

Elastic interaction forces are proportional to the disk 
overlap, /fj — k{Ri + Rj — rij)"-<'\ where fij := fj — fi 
points from one bubble center to another and Ri labels 
the radius of disk i. In the full model that we focus on 
here, viscous forces oppose the bubbles' relative veloc- 
ity Avij := Vj — Vi with magnitude f^-^'^ = 6|At;ij|"^''"^ 
|12| . In simulations we set both exponents ad and ctvisc 
to unity. The strain rate 7 is imposed via Lees-Edwards 
boundary conditions. The unit cell contains a 50:50 bidis- 
perse mixture of TV = 1020 to 1210 bubbles with size ratio 
1.4:1. Stresses are averaged over a run (total time 2O/7) 
after discarding the transient. 

Elastic and Viscous Stress — Because forces balance 
on each bubble, the shear stress can be computed accord- 
ing to a^y := CTtot = 217 E(ii> rtj,x (ffly + ftj^y), where 
V is the area of the unit cell and the sum runs over 
contacting pairs. It is convenient to distinguish contri- 
butions of elastic and viscous forces to the total stress 
CTtot := o-ei -I- (Tvisc- Fig. [T] dcpicts aei and a^isc as func- 
tions of strain rate for three packing fractions. We find 
that the viscous stress Cvisc depends only weakly on 
and scales linearly with 7, and dominates the total stress 
for strain rates 7 > 0(10^^). We denote this regime as 
Viscous (V), and since (ivisc ~ 7, crtot ~ 7 here. Our ana- 
lytical model, developed below, treats the case where CTci 
dominates the stress, and in all that follows, 7 < 10~^, 
CTvisc ^ CTei, dtot ~ CTci, and we no longer distinguish 
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FIG. 2: (a) The probability distribution function (pdf) of 
relative bubble velocities |At;|, for <^ « (/)c and strain rates 7 
spanning three decades (see legend), does not collapse when 
rescaled by 7. In contrast, the pdf of |A«|/7^'''* yields a rea- 
sonable collapse (inset). Note that Eq. ([TJ predicts the pdf's 
second moment, not its shape, (b) Collapsed stress-effective 
strain relation for 7cff oc -f 0.7(7/0-)^''^ and the same data 
as in Fig. [i] The sohd curve is j/ = O.OSSxVl + O.OSa;^. 

between utot and (Td, referring to both as as a. 

Main Phenomenology — As Fig.[l|D illustrates, the rhe- 
ology is nontrivial and departs from simple linear scal- 
ings. For above (^^ = 0.8423 ± 0.001 [I3], the stress 
flattens as 7 is lowered, while for « 0^ the stress is well- 
described by a power law g ^ 7^/2 — this will emerge 
below as the Critical scaling regime of the rheology. For 
(j) below (f)c, the stress shows increasing downward curva- 
ture as strain rate is decreased; our model does not treat 
this case, and we do not consider it further. 

Analytical Model — To relate stress and strain rate, 
we will (i) find expressions for the dynamic and yield 
strains 7y and 7dyn that constitute the effective strain 
and connect them to the stress via power balance; (ii) 
express the stress as a function of 7eff . 

(i) The effective strain jdyn can be thought of as the 
typical strain undergone between plastic rearrangements 
of the contact network. Because bubbles move with rel- 
ative velocity |Au|, the contact network rearranges on a 
time scale ^dyn ~ c^/IAul, where d is the average bubble 
diameter. The typical strain incurred on this time scale 

is 7dyn := 7tdyn ~ 7/|Al>|. 

To obtain an estimate for |A?;|, we turn our attention 
to the relation between the scaling of stress and velocity 
fluctuations. As noted in Ref. [M], mechanical energy is 
supplied to the system at a rate ~ crtot7- Energy dis- 
sipation takes place by bubbles moving past each other 
— hence the dissipation rate scales as /^'^lAul ~ |Aup 
|15j . Balancing the two yields 
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Eq. ([T]) is the first of three relations comprising our 
model: given the stress, it provides the scaling of |Aw|. 

It will emerge from our model that the nontrivial scal- 
ing of velocity fluctuations underlies the rich rheology. 
Fig. [2^ shows probability distributions of IAwI/t- for 
(j) ^ (j)c and 7 < 10~^. Since here tr ~ 7^^^, |Au| does 
not scale as 7. In fact, Eq. ([T]) predicts that in this case 




FIG. 3: Schematic depiction of predictions for scaling of stress 
a and fluctuations in the relative bubble motion |Aw|/7 with 
strain rate 7 and distance to the critical packing fraction 
Atp = (j> — (pc- There are four distinct flow regimes: Yield 
Stress (YS), Transition (T), Critical (C) and Viscous (V). 

I Aw I ^ j-^/^. This gives a good collapse of the data (inset 
Fig. [2^). Note that in the Viscous regime where cr ~ 7, 
one finds that \Av\'^ ~ 7^, so that the typical relative 
velocity |Au| scales trivially with 7 (not shown). 

For 4) > (t)c one anticipates a threshold (yield) stress 
even for vanishingly low strain rate, which in our picture 
translates into an additional contribution to the effective 
strain; this is 7y. A reasonable expectation for the scaling 
of 7y is the strain scale required to prepare a packing at 
4> = 4'c ~\~ A(^ by compressing a system from the critical 
packing fraction: 7y ~ A0/0 ~ Collecting terms, 

the effective strain reads: 



7eff = yliA0 + A2d7/|Aw| 



(2) 



(ii) We now construct a stress-strain relation a = 
g{A.(j)^ 7cff ) 7cft and make the ansatz that the shear mod- 
ulus g displays single parameter scaling: g{lS.4>,^cs) — 
A0P g(7off/A0'). We will determine a form of g based 
on known results for static systems. Above (pc, static 
systems display a regime of linear response a = Gj, 
where the static shear modulus G = Goy/Acf) [TB]. Hence 
p = 1/2 and g{x) — > Go for x 0. Precisely at (j)c, G 
vanishes, no analytic expansion of the stress-strain rela- 
tion is possible, and critical static spring networks dis- 
play the quadratic form a = K|7cff|7eff [IZj- It follows 
that g = 1/2 and g{x) — > k\x\ for x — 00. Therefore the 
stress-strain relation can be rewritten as 
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In the analysis to follow, only p, q, and the asymptotic 
scaling of g{x) are essential. In Fig. [2]d we show that 
a scatter plot of ct/A0 as function of 7eff / V A0 shows 
excellent data collapse with the correct asymptotic be- 
havior — here jcS oc A(l>+{A2d/Ai){'j/a)'^/'^ and A2/A1 
has been adjusted to obtain collapse. We also note that 
the simplest choice for g{x) that obeys reflection symme- 
try, remains analytic above (jic and obeys all necessary 
Gov^ 



scalings is g{x) 
remarkably well (Fig. [2|3). 



Kx/Go)^, which fits the data 
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FIG. 4: Scaling collapse over four decades in strain rate 7 
and three decades in distance to the critical packing fraction 
A(j> (legend). Rescaled coordinates a/A(f> and 'y/A(f>'^ are ap- 
propriate for parameters spanning the Transition and Critical 
regimes (see Table Dashed lines are guides to the eye with 
slopes 1/3 and 1/2. (inset) Boundaries between the Yield 
Stress, Transition and Critical regimes in the A<j)--y plane. 



TABLE I: Analytical model and its solutions in the Yield Stress, Transition, Critical and Viscous regimes. 

This strong test of the model is shown in Fig. |4j where 
data for three decades in A0 and four decades in 7 col- 
lapse to a single master curve. Data points satisfying 
7 < ciA^'''/^ and 7 > [18i are labeled Yield Stress 

(black) and Critical (blue), consistent with scaling pre- 
dictions; the Transition regime (red) lies in between. YS 
data points "peel off" from the master curve (note the 
black data points above red ones) because it is not pos- 
sible to collapse three regimes in one plot when their 
crossovers scale differently. 

We stress that the scaling exponents, including depen- 
dence on A(f), are all predictions, not adjustable param- 
eters. The excellent data collapse in Figs. [2]d and |4] is 
therefore a striking confirmation of the model. 

We can gain some intuition for the various regimes by 
considering different approaches to the critical point (see 
Fig. [4] inset). Fixing (f> = (j)c and adiabatically lower- 
ing the strain rate approaches point J from the Critical 
regime a ^ 7dyn ^ 7^^^i where stress is always domi- 
nated by dynamic effects. Similarly, fixing 7 = 0"*" and 
adiabatically decreasing Ac/) approaches point J from the 
Yield Stress regime a ^ Gjy ^ A(/)^^^, where the flow 
is rate independent. Finally, there is an anomalous flow 
regime a ~ G'7dyn ~ A(/)^/^7^/'^ that transitions between 
the Critical and Yield Stress regimes. It is traversed when 
varying 7 at finite A4> or vice versa. 

Length scale — In our model there is no strain rate 
dependent length scale [H [TU] to capture "swirls" or 
avalanching rearrangements. To test this assumption, in 
Fig. [5^ we plot correlations in the non-afhne component 
of the bubble velocities, C{y) := {vx{0)vx{y)) ■ We find 
C{y) to have a form that is (i) reminiscent of disordered 
static [12] and quasistatic [7] soft sphere systems and {ii) 
remarkably robust to changes in 7 and (j). Moreover, as in 
static linear response [12], C{y/L) collapses for different 
box sizes L, suggesting that for the system sizes studied 
here the box size is the only relevant macroscopic length 
scale. We note that C{y) was measured in Ref. [1] in the 
mean field (MF) bubble model, which replaces realistic 
bubble-bubble viscous forces by an effective drag term 
that punishes deviations from an afhne (linear) velocity 
profile — see Fig. [sja. The minimum of C{y) in the MF 
model (Fig. [5^ inset) selects a length that was found to 
scale with A0 and 7 Because this behavior vanishes 
when the mean field approximation is lifted, we conclude 



Our model comprises Eqs. (T][3|, which express a, jcS, 
and |Ati| in terms of 7 and Acj). For scaling analysis the 
constants Ai^ A2, Gq, k and d can be set to unity. 

Flow Regimes — The three equations for a, 7off, and 
I Aw I in terms of 7 and Acj) lead to our scaling predic- 
tions. Eqs. ([2]) and ([s]) each have two scaling regimes, 
which are selected by varying 7 and A(j). In combina- 
tion, these contribute three scaling regimes to the rheol- 
ogy, i.e. the stress-strain rate relation. So where previous 
scaling ansatze presume two rheological regimes [4j |6l |8] 
— a yield stress plateau giving way to a power law in 7 for 
higher strain rates — we find Yield Stress (YS), Transi- 
tion (T), Critical (C) and Viscous (V) regimes, each per- 
sisting over a finite range of strain rates (Fig.jsj). Table I 
collects the pertinent equations, solutions and parameter 
ranges for all scaling regimes. 

Checking the Model — For data restricted to two 
regimes, it is possible to collapse the j-a flow curves to 
a master curve by rescaling with At/). For the Transition 
and Critical regimes, the rheology is predicted to obey 
a A(/)3 73 and ct ^ 72 , respectively. Hence for data in 
these two regimes, a/A(f> vs j/Acfy^ can be collapsed to 
a master flow curve, characterized by a crossover from a 
1/3 to a 1/2 power law scaling. 
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FIG. 5; (a) Two-point correlation function C{y) for L — 75 
(as in Figs, [l] and [i]), 7 = 10'^ and 10~^ and = 0.82, 
(j> ~ (jic and = 0.86. Also plotted are system sizes L = 54 
and 105 for strain rates 7 = 10~^ and 10"* and packing 
fractions (j) = 0.82, 0.84, and 0.86. (inset) C{y) has a different 
form in the mean field model and evolves with 7 and A(^. 
Here, two examples: </> = 0.841 (solid) and 0.88 (dashed) at 
7 = 10~^. (b) Dissipation in the models. Viscous forces 
(black arrows) are proportional to the difference between a 
bubble's own velocity and a linear velocity profile (mean field 
model) or the velocity of a second bubble (full model). 

that it is an artifact of the MF dynamics. 

Discussion and Outlook — The unanticipated presence 
of three rate dependent regimes with distinct rational- 
valued scaling exponents offers an explanation for dissim- 
ilar exponents found in the literature [iMl IFl [TUl [TTl ED] . 
For data spanning some combination of the Transition, 
Critical and Viscous regimes, one could fit a power law 
cr ^ 7^ with effective exponent 1/3 < (i < 1. Re- 
ported scalings indeed range from around f3 « 0.4 [4 
to 0.6 0(20] and even 1.0 HI]- Though the prediction 
is more difficult to test, the model also provides a plau- 
sible argument that the dynamic yield stress of systems 
with spring-like elastic interactions has super-linear de- 
pendence on distance to (j)c (cy ^ A0^/^), rather than 
linear [HI HI]- Refs. [HH] indeed find super-linear scaling, 
while (Ty A(j) in the quasistatic simulations of Ref. [7] . 
Recent work by Hatano probing the lowest shear rates to 
date finds ay - Ac/)^ with A = 1.5 ± 0.1 pSl. 

The model is easily generalized to other microscopic in- 
teractions. Elastic interactions enter through the scaling 
of the shear modulus G ^ A^""'^^/^, while different vis- 
cous force laws affect the fluctuations via power balance: 
CT7 ^ I Au|"^'=<^+^. Recent data for viscous NIPA particles 
are consistent both with our prediction of cr 7^/^ in the 
Critical regime and with a yield stress dy ~ A^^^'+^Z^ for 
Hertzian interactions, ad — 3/2 [22 . For physical foams, 
beheved to have a viscous exponent ctvisc = 2/3 (Ml, the 
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Critical regime scales as cr ~ 7 
in remarkable agreement with recent experiments that 
found a ~ -jO-se Finally, for slow frictional flows, 

both the the drag forces and the global rheology are rate 
independent [3S]. We suggest that the global rate inde- 
pendence is not a triviality and note that it is consistent 
with our model, where cr ~ 7° for ayisc — ^ 0. 
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